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Abstract 

The photoluminescence decay in amorphous semiconductors is described by power law t~ s at long 
times. The power-law decay of photoluminescence at long times is commonly observed but recent 
experiments have revealed that the exponent, 5 ~ 1.2 — 1.3, is smaller than the value 1.5 predicted 
from a geminate recombination model assuming normal diffusion. Transient currents observed 
in the time-of-flight experiments are highly dispersive characterized by the disorder parameter a 
smaller than 1. Geminate recombination rate should be influenced by the dispersive transport 
of charge carriers. In this paper we derive the simple relation, 5=1 + a/2. Not only the 
exponent but also the amplitude of the decay calculated in this study is consistent with measured 
photoluminescence in a-Si:H. 
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I. INTRODUCTION 



The photoluminescence in amorphous semiconductors consists of a major band and 
weaker bands of longer wavelength attributed to dangling bonds jl| . The main photolumi- 
nescence band arises from the radiative recombination of electrons and holes both trapped 
in band tail states Our interest is in this main band excited with relatively weak light 
pulses at the absorption edge. Photoluminescence decay is described by power law t~ s at 
long times above a certain temperature . For sufficiently law excitation 

intensity a geminate recombination model for photoluminescence decay has been developed, 
which predicts a t~ 3 ^ 2 decay of luminescence at long times In the model the t ' 

behavior at long times is attributed to normal diffusion of a geminate pair. Predicted t~ 3 / 2 
decay of luminescence at long times is confirmed by early experiments, where high tempera- 
ture data is different from the low temperature decay for which diffusion is negligible y], 0] ■ 
At high temperatures electrons have enough thermal energy to conduct hopping transport 
in the band tails l], 0, 0] . It is widely recognized that a hole is trapped at a deep localized 
very short times after the pulse excitation The power-law decay of photolumi- 
nescence at long times is commonly observed but recent experiments have revealed that the 
exponent is smaller than the predicted value, 5 ~ 1.2 — 1.3 < 1.5 □ □□□□ . The exponent 
increases with increasing temperature. In other words, the deviation of the exponent from 
1.5 decreases by raising temperature. Such deviation is interpreted in terms of deviation of 
motion of charges from the normal diffusion DDDOD 

Indeed, at temperatures at which power law of the photoluminescence decay is observed 
with 5 smaller than 1.5 transient currents observed in the time-of-flight experiments are 
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B hly disp ersive and completely different from those obtained from normal diffusion 
ll2L Il3l. The shape of the transit current, /, is featureless in the conventional units / 
versus t 10, The transit time, t tr , is obtained only from the log/ versus logt 

traces, because many time scales coexist for electron hopping between localized sites in the 
band tail states jll|. The shape of the transit current is derived by a continuous time random 
walk model where the hopping time distribution function has an algebraic asymptotic form 
jiol 11 . 14, 15| ip(t) ~ ^- Q £-( 1+Q ) ; where 7 r is a carrier jump frequency and a is a constant 



in a range of < a < 1. From the current traces the disorder parameter, a, is typically 
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obtained as a ~ 0.6 0, HQ, 0, Q. It implies that the displacement, r, of a charge 

carrier obeys (r 2 (t)) oc t a . Therefore, electron transport cannot be described by normal 
diffusion for which (r 2 (t)) oc t. 

Geminate recombination rate should be influenced by dispersive transport of charge car- 



riers 



The fact is widely recognized but a simple relation between 5 and a is 
not established. In this paper we derive a simple relation between the long time exponent 
of photoluminescence decay, S, and the disorder parameter, a. The relation is confirmed by 
experimental data of a-Si:H taken at various temperatures. In the experiments, the samples 
of a-Si:H were prepared by a glow discharge method. The sample held in helium gas was 
excited with 10 nsec light pulses of energy around 2.0 eV at the absorption edge in a-Si:H. 
The photoluminescence main band is observed at the peak energy 1.22 — 1.35 eV. We also 
obtain a theoretical relation between amplitude of the power law decay of photoluminescence 
and the disorder parameter, a. The amplitude decreases due to escape of an electron-hole 
pair from recombination. The result is also confirmed by measured photoluminescence data. 

II. PHOTOLUMINESCENCE 

Amorphous semiconductors have localized states in band tails jjj. The localized states 
are randomly distributed and their electrostatic energies are also distributed. In a-Si:H 



y trapped and electrons 
l|. Radiative recombina- 



electron-hole pairs are created by light irradiation. Holes are quic 
execute hopping-random walks among localized band tail states 
tion of geminate electron-hole pairs is observed as photoluminescence whose decay curve 
is independent of the excitation intensity in a range from 5[nJ/cm 2 ] to 500[//J/chi 2 ] jsj]. 
Intensity of photoluminescence is proportional to the recombination rate which is given by 
the decay rate of survival probability of electrons. Long time decay of photoluminescence is 
controlled by diffusion of electrons to the sites where holes are located. Motion of electrons 
due to hopping-random walks has been measured by transient current in the time-of-flight 
experiment |l0j, 
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displacement of electrons obeys 

(r(t) 2 )<xt a . (2.1] 
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The initial distance of electron-hole pairs may have spatial distribution even immediately 
after photo-excitation. However, the long time asymptotic kinetics of photoluminescence 
decay is insensitive to the initial distance, r . The distribution of initial distance immediately 
after photo-excitation is important to analyze short time kinetics of photoluminescence. In 
this paper we focus on the long time asymptotic kinetics. As written before, intensity of 
photoluminescence is proportional to the recombination rate. In the limit of infinitely fast 
intrinsic recombination rate (recombination rate at an encounter distance) compared to the 
thermal jump frequency of electrons the recombination reaction rate, Rcx>{^o,t), is nothing 
but the probability distribution of first arrival at the encounter distance R; namely, the 
first-passage time density. The transition probability or the conditional probability g (r ,t) 
of finding an electron at the encounter distance at time t when it starts initially at site r 
in the system with no reaction satisfies the integral equation 17 



g{f ,t)= [ dti0(iM-ti)iMn>,ti) (2-2) 
J o 

because random walker at the encounter distance at time t must have been there for the 
first time at time t\ < t. After Laplace transformation of Eq. (|2.2|) Roo (r , s) is obtained as 

In experiments the intrinsic recombination rate is not infinitely large but it competes with 
jump frequency of electrons. By increasing temperature the luminescence quantum efficiency 
decreases, which indicates that more electrons escape recombination due to increase of hop- 
ping frequency. The effect is strong when the hopping frequency competes with the intrinsic 
recombination rate. When recombination rate competes with jump frequency of electrons, 
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recombination rate, R (r , t), is obtained from R^ (r , t) by the method of Pedersen 

First we note that the probability that a particle which starts at r will visit the spherical 
shell of radius R during the time interval between t and t + dt is equal to the recombination 
rate for the perfectly absorbing boundary condition. Two types of waiting time distribution 
is defined inside the reaction zone; one is the waiting time distribution function of making 
a jump without reaction, ip ut{t), and the other is the waiting time distribution function of 
reaction, ip rc {t). The waiting time distribution of making a jump at the encounter distance 
is given in terms of that in the absence of reaction, ip(t), as, 

V>out(t) = exp (- 7rc t) . (2.4) 



The waiting time distribution of reaction is denned as the probability that the particle which 
is initially at a site in the encounter distance will undergo reaction without making a jump at 
time t. It is given by the reaction rate constant, 7 rc , multiplied by the remaining probability 
of particles at the site in the reaction zone, which decays either by jump motion or reaction, 



iprc(t) = 7r C exp (-7 rc t) / dt^fa). 



(2.5) 

Some particles may recombine at the first visit to the encounter distance R. Other fraction 
of particles escape recombination at the first encounter and make a jump to the spherical 
sphere of radius R + b with b being the jump length. Some particles may recombine at the 
second encounter after jumping from the sphere of radius R + b. The recombination rate 
R(fo,t) should be the sum of the contributions from the first encounter, second encounter 
and so on, 



R{r ,t) 



t 

dhtprc (t - ti) Roc (r ,ti) + 



dt x dt 2 I dt 3 iP IC {t-t 1 )R oo {R+b,t 1 -t 2 )iP mt {t 2 -t 3 )R oo {f ,t 3 ) + --{2.6) 
o Jo Jo 



By Laplace transformation Eq. ()2.6|) becomes 21 1 

R(r ,s) = R, 



\r ,s) 



4>r 



(2.7) 



1 - tpoMt(s)Roc (R + b,s) 

Therefore, R (r*o, s) is obtained from R^ (r*o, s), which is expressed in terms of the transition 
probability g (r , s) as Eq. f)2.3j) . When migration of charges is due to thermal transition be- 
tween trap states and is described by fractional Brownian motion, the transition probability 



is given by 
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exp 



9 (r , s) 



r -R 



D„s 



4irD a roS 



l-a 



R 



(2.8) 



where the perfectly reflecting boundary condition is imposed at r = R. When the fractional 
Brownian motion is derived from the continuous time random walk with jump length b and 
jump frequency 7 r , a generalized diffusion coefficient is derived 0,13; 



D, 



sin ira 
2na 



ai2 



7 r Q 6 



(2.9) 



from the precise form of the waiting time distribution function of the release from the trap in 

nn 

the absence of recombination reaction IU,115|], which is written after Laplace transformation, 

fa)~l-J^(±Y. (2.10) 

SUl 7TQ! \7r/ 

By substituting Eq. ()2.8|) into Eq. ()2.3|) we get i?^ (r ,s), which is related with R(r ,s) 
through Eq. ()2.7|) . The recombination rate in Laplace domain is derived as, 



R (f , s) = — exp 
r 



r -R 

jD n s~ c 




R 



(2.11) 



V D a s~ 

where a generalized intrinsic reaction rate is defined as |24j, |2| 

a sin ircy 

k a = 4j vc {0)- -^2nR 2 b ~ ^ c 27rR 2 b. (2.12) 

na 

After inverse Laplace transformation the asymptotic photoluminescence decay is obtained 
as 

R(- t \ a R ( r *l R ~ l , 4*RDa/k a \ R 1 , , 

2T{l-a/2)r»\l + ^RD a /k a + (1 + A-nRDjKf) v^t 1+a/2 ' K } 

According to Eq. ()2.13|) the power of asymptotic photoluminescence decay, l/t s , is related 
to the disorder parameter a by 

5 = | + 1, (2.14) 

for any strength of the intrinsic recombination rate, k a , and the initial distance between the 
electron and the hole of a geminate pair. In Fig. 1 the exponent 5 of the long time power-law 
decay of photoluminescence is shown together with the exponent 1 + a/2 calculated from 
the disorder parameter, a, of the dispersive transient current, which is measured by electron 
time-of-flight experiments at low electric fields [l^ . The two exponents are consistent, which 
indicates that the long time photoluminescence is controlled by the dispersive transport of 
electrons. 

The intensity of photoluminescence depends on the strength of intrinsic recombination 
rate and the jump frequency of electrons. The intensity of photoluminescence decreases 
by increasing temperature because the number of electrons which escape recombination 
increases. We neglect other non-radiative processes which are important in samples with 
a substantial defect density [l]. Obviously, the intensity is also influenced by the initial 
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distance between the electron and the hole of a geminate pair and the intrinsic recombination 
rate. We investigate two extremes, separately. The first one corresponds to the band- 
edge excitation, where an electron is created in the vicinity of a hole and the intrinsic 
recombination rate competes with the dispersive transport of electrons. The second one 
corresponds to excitation of electrons to the conduction band through which an electron is 
separated from a hole by large distance and the recombination occurs whenever the electron 
encounters the hole. 

In the first case the initial distance between an electron and a hole is small and intrinsic 
recombination rate k a is small. An electron may visit the encounter distance many times 
before recombination. In experiments the wavelength of incident light is chosen to excite 
electrons in the tail of band edge and electrons are excited close to the immobile holes. For 
simplicity the initial distance of an electron-hole pair is assumed to be in the vicinity of the 
encounter distance, R (R + e, e — > 0) . The recombination rate is obtained from the inverse 



Laplace transform of Eq. (j2.11)) as 



R{r ° = R ' t) = 4^DJ E *M 2 I' AnRD a /k a R~ ' ' ^ 



where the generalized Mittag-Leffler function is defined as 

OO 

^ w= Sw (2 ' 16) 

In the long time limit the above equation reduces to 

3/2 



R ^ = ^) ~ f (j^) r(l-a/2) 



3/2 I ^l+a/2 



(2.17) 



(7 r 3/ 77rc)V +a/2 ' 

where C\ is a constant. The intensity of photoluminescence decay depends on the relative 
magnitudes of 7 r and 7 rc . The photoluminescence is weak when the jump frequency of 
electron migration is large because electrons escape recombination by diffusion. Strong 
photoluminescence is expected when the intrinsic recombination rate is large. 

In the second case the initial distance between the electron and the hole of a geminate 
pair is assumed to be large and the intrinsic recombination rate is assumed to be infinitely 
large. Then we obtain from Eq. (|2.13|h 



where Ci is a constant 27[. Intensity becomes smaller as the jump frequency increases. 



In both cases the photoluminescence intensity has a scaling form, 

(2-19) 

where C is a constant. The photoluminescence decay measured by Murayama et al. obeys 
the scaling given by Eq. (|2.19|) . if temperature dependence arises only from the disorder 
parameter, a, and values of a measured by the time-of- flight experiments at each tempera- 
ture are introduced. The temperature dependence of a is phenomenologically described by 
a linear function, 

a = a + ^, (2.20) 

as shown in Fig. 1, where = 0.15 and T = 313.4 [K]. If we assume an exponential 
distribution of band tail states characterized by an attenuation parameter k^T , a is given 
by a = T/T . The presence of the constant term, ao, in Eq. (J2.20|) implies that activated 
release of an electron from an exponentially distributed band tail states is a too simplified 
model. Nevertheless, the linear function is a good approximation for observed values of the 
disorder parameter, a. By introducing the linear temperature dependence of a, Eq. ()2.20|) . 
into the scaling relation, Eq. (|2.19|) . photoluminescence intensity at a certain time, t = to, 
after photo-excitation is rewritten in the following way as a function of temperature, 

R(f ,t )~exp(-pj, (2.21) 

where T c = Tq/ In (Cy/to). In Fig. 2, photoluminescence observed at a delay time of to = 10 
[//sec] is plotted against temperature. The line represents Eq. ()2.21|) with T c = 30 [K]. 
The similar temperature dependence is observed for quantum yield of photoluminescence 
1,I3|. Since the kinetics of photo-luminescence decay is highly non-exponential, quantum 
yield does not necessarily obey the same temperature dependence as photoluminescence 
intensity at a certain delay time. Eq. (|2.21|) is proposed for quantum yield by assu ming 
activated behavior for the competition between radiative and non-radiative processes [29]. 
On the other hand, in our derivation of photoluminescence intensity at a certain delay 
time, escape of an electron from recombination is a non-radiative process which competes 
with radiative recombination. By substituting experimental values of T , T c , and t into 
the relation, Cy/to = exp(To/T c ), we find, C 2 = 6.8 x 10 13 [l/s]. Now, we consider two 
extreme situations leading to the same scaling form Eq. (|2.19|) . Although the temperature 
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dependence is the same, different values of parameters may be assigned in these two cases. 
In the first case we find 7 r ~ 10 14 [l/s] for 7 r ~ 7 rc [l/s]. The values are consistent with 
7 r ~ 10 14 — 10 15 [l/s] estimated from the relation {L/b) 2 ~ (7 r ^tr) Q ; where a ~ 0.6, the zero- 
field value of the transit time, t tT ~ 10 _4 [s], in time-of-flight experiments with the sample 
thickness, L ~ 5[/xm], and the jump length, b ~ 20 — 50A, assumed for hopping random walk 
of an electron are introduced. In the second case we obtain 7 r ~ 6.8 x 10 13 [l/s] which is 
also consistent with the value estimated above. The measured photoluminescence obeys the 
scaling form, Eq. (|2.19|) . with reasonable values of parameters for each case. Unfortunately, 
it is impossible to differentiate the results of these two extreme situations only from the 
scaling relation of long time asymptotic. 

In order to see the time range over which the long time asymptotic, Eq. (J2.19|) . is valid, 
lines given by Eq. (j2.19|) are compared in Fig. 3 with measured photoluminescence decay. 
It is difficult to keep all experimental conditions unaltered when temperature is changed. 
Therefore, intensity difference as a function of temperature is not reliable and we focus our 
attention on its time dependence; we used Eq. (j2.19j) with temperature dependence of a 
given by Eq. (J2.20j) and C as a fitting parameter. As we can see from Fig. 3, the lines 
given by Eq. (j2.19|) are consistent with experimental results at long times at relatively high 
temperature, T > 105 [K]. However, short time kinetics deviates from long time algebraic 
decay. The transition occurs at larger values of time for lower temperature. Below 100 [K], 
the photoluminescence decay at to = 10 [//sec] is not in the regime of long time asymptotic. 
This is the reason for the deviation from the scaling relation at temperatures below 100 [K] 
in Fig. 2. At low temperatures electron-hole recombination tends to occur before electron 
executes random walk by thermal hopping. Such short time kinetics is out of the range of 
applicability of the present model and requires extending our theory to include the precise 
distance dependence of the intrinsic rate. Nevertheless, the observed slow decay kinetics 
up to 10 microseconds at 105 [K] in Fig. 3 suggests that the intrinsic reaction rate has a 
finite value and competes with hopping transition. Thus, we can draw the following picture. 
Initially an electron is photo-excited close to a hole by band edge excitation. The hole is 
immediately trapped. The electron executes random walk by thermal hopping. It is most 
likely that the electron does not necessarily recombine with the hole even when it reaches a 
recombination distance 28^ . 

The recombination rate in a lattice model is also derived, where an electron is excited 
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initially at a neighboring site of a hole and recombination takes place with a certain rate 
when an electron and a hole occupy the same site. The recombination rate is obtained from 
the decay rate of survival probability, N lbj,t), calculated previously as, 



N[h " t] « WiW^ ' (2 ' 22) 



which is Eq. (6.13) of ref. 2l|. By taking the time derivative of the above equation the 



scaling relation, Eq. (|2.17|) . is reproduced. 



III. DISCUSSION 



The photoluminescence observed at temperatures higher than 50 K has a long time 
asymptotic decay described by t~ 5 . The exponent, 5, is smaller than 1.5 predicted from 
a geminate recombination model with normal diffusion of an electron If one takes 
into account that electron hopping between localized sites in the band tail states cannot 
be described by normal diffusion but by dispersive transport due to many time scales for 
detrapping, the exponent 5 is expressed in terms of disorder parameter, a by 5 = 1 + a/2. 
The predicted value of S is consistent with experimentally measured value when a obtained 
by time-of-flight experiments is used. In the time-of- flight experiments the transit current, 
/, is described by two different power law decay functions, / oc t~( 1_Qi ) (£ < t tr ) and 
/ oc t~( 1+Q ^') (t > t tr ), where t tr is called the transit time. The continuous time random 

nn 

walk model and fractional Fokker-Planck equation predict a« = aj [10, I27| . However, 
measured values of a, and aj are slightly different 3(|. When they are different we 
choose a = oti because photoluminescence mainly decays within the transit time and the 
origin of discrepancy between ai and a/ has not been fully understood. A possible origin is 



the structured distribution of localized states 
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32]. 



The other possibility is that it is due 



to the correlation between waiting times which is neglected in the continuous time random 
walk model. In fact different values of on and aj are obtained in the hopping simulation in 
which the correlation between hopping frequencies exists .3.3] ■ 

Theory and experiment can also be compared by evaluating the temperature dependent 
decrease of the luminescence intensity due to escape of an electron-hole pair from recombi- 
nation. The measured intensity of photoluminescence obeys a scaling relation, Eq. ()2.19|) . 
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The scaling relation is obtained from a geminate recombination model with dispersive dif- 
fusion. The reasonable parameter values are obtained by assuming that initially an electron 
is photo-excited close to a hole and the intrinsic recombination rate of an electron with a 
hole is not infinitely large. 

The competition between radiative recombination and non-radiative escape processes of 
diffusion is taken into account in geminate recombination model by introducing a sink term 
for radiative recombination and dispersive diffusion for transport. Recently, we have derived 



a fractional reaction-diffusion equation from a continuous time random walk model 
Recombination rate, Eq. ()2.13j) . can be also derived from the fractional reaction-diffusion 
equation, 



d . s d /•* 1 1 



D a V 2 p (r, h) - k a K ^ R2 > p (r, h) 



(3.1) 



dt ry ' ; dtj x r(«) (t -to 1 

Another remark on the temperature dependence of the photoluminescence intensity is that 
measured values of a is not proportional to temperature for some semiconductors |3^| . In 
principle, the measured values of a for each temperature should be employed when we 
compare theoretical results with experimental data. 

As in the previous geminate recombination model for normal diffusion, we have neglected 
the coulombic interaction of an electron with a hole, althoug 
phous semiconductors are known to be large close to 100A 



i Onsager radius in these amor- 
For normal diffusion, the 



power of 



ong time algebraic decay of recombination rate is independent of coulombic inter- 



action . In the presence of distributed site energies of localized states, we should further 
take into account the relative magnitude of the coulombic interaction against the fluctua- 
tion of activation energy associated with hopping transition. When an electron is distance 
r away from a hole coulombic potential is U(r) = — e 2 /(er), where e ~ 10 is the dielectric 
constant of a-Si:H P|. Coulombic interaction creates energy difference for a hopping transi- 
tion with a distance b ~ 20 — 50A, which is estimated to be bdU(r) /(dr) = e 2 b/(er 2 ), while 
the fluctuation of activation energy associated with a hopping transition is characterized by 
the mean depth of distributed site energies of localized states, k-gTo with Tq ~ 313 [K]. By 



equating these two quantities, we find a characteristic length, R s = ^/e 2 b/ (eK-bTq) ~ 40A, 
which is much smaller than the Onsager length. When distance between a hole and an 
electron is larger than the characteristic length, R s , coulombic interaction can be neglected. 
Since R s is small and comparable to the hopping distance and we are interested in the long 
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time asymptotic of photoluminescence, where electron migration at large distances from a 
hole is important, the power of the long time asymptotic of recombination rate would not 
be influenced by coulombic interaction. The interaction may be important for the decay of 
photoluminescence at short times, which we do not investigate in this paper. Further theo- 
retical study is required to clarify the influence of coulombic interaction on recombination 
rate. 

A probable mechanism of geminate recombination is radiative tunnelling which occurs 
at short distance comparable to hopping distance . Such short range reaction can be 
modeled by recombination reaction at encounter distance as long as long time asymptotic 
is concerned. Distance dependence of intrinsic reaction rate matters for short time kinetics. 
The more precise model should be introduced to account for the photoluminescence decay 
at short times. Murayama et al. has developed a theoretical model for dispersive tr ansp ort, 
where the energy fluctuation of the localized band tail states is self-affine 7, 16,0, [JS]- 
By the Monte Carlo simulation the photoluminescence decay over the whole range of time 
are well reproduced as well as the dispersive transit currents. More recently, the fractional 
diffusion equation is numerically solved by Fukunaga to reproduce the photoluminescence 



decay over the whole range of time 



although his way of including recombination reac- 



tion is different from ours because the reaction term is not coupled to the memory kernel like 

Eq. dsnj 



25l l3q . The long time asymptotic is different from ours when a > 2/3 21 1. 
Although our theory is restricted only to the long time asymptotic decay, it provides an ana- 
lytical relation between the exponent of algebraic photoluminescence decay and the disorder 
parameter a, which can be obtained from the time-of- flight experiments. The amplitude of 
the power law decay of photoluminescence is also given in terms of the disorder parameter, 
a. Some new experimental evidence is also presented in this paper, which supports our 
analytical results. 
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Figure 1: The exponent 5 of the long time photoluminescence decay and the disorder 
parameter a, in a-Si:H against temperature. Closed triangles are 5 obtained from the total 
light decay of photoluminescence in a-Si:H excited with the pulsed laser of photon energy 
of 2.32 [eV]. Closed circles are a obtained from electron time-of-flight experiments in a-Si:H 
jl6^ . Open circles are calculated by 5 = 1 + a/2 from values of a shown by closed circles. 
The solid line is a linear fit, a = 0.15 + T/313.4. 

Figure 2: Temperature dependence of photoluminescence intensity from a-Si:H excited 
with the pulsed laser of photon energy of 2.32 [eV]. The photoluminescence has been observed 
at the delay time 10 [/xsec] and the peak energy of 1.24 [eV]. Solid lines is exp (— T/T c ) with 
T c = 30 [K]. 

Figure 3: Total light decay of photoluminescence in a-Si:H excited with the pulsed laser 
of photon energy of 2.32 [eV]. The photoluminescence has been observed at different tem- 
peratures, 105 [K] , 145 [K] and 170 [K] (top to bottom) j^j]. Solid lines are theoretical lines 
given by Eq. (12~THD . 
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